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Exploring star formation using the filaments in the Sloan 
Digital Sky Survey Data Release Five (SDSS DR5) 
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ABSTRACT 

We have quantified the average filamentarity of the galaxy distribution in seven nearly 
two dimensional strips from the SDSS DR5 using a volume limited sample in the ab- 
solute magnitude range —21 < Mr < —20. The average filamentarity of star forming 
(SF) galaxies, which are predominantly blue, is found to be more than that of other 
galaxies which are predominantly red. This difference is possibly an outcome of the fact 
that blue galaxies have a more filamentary distribution. Comparing the SF galaxies 
with only the blue other galaxies, we find that the two show nearly equal filamentarity. 
Separately analyzing the galaxies with high star formation rates (SFR) and low SFR, 
we find that the latter has a more filamentary distribution. We interpret this in terms 
of two effects (1.) A correlation between the SFR and individual galaxy properties 
like luminosity with the high SFR galaxies being more luminous (2.) A relation be- 
tween the SFR and environmental effects like the density with the high SFR galaxies 
preferentially occurring in high density regions. These two effects are possibly not 
independent and are operating simultaneously. We do not find any difference in the 
filamentarity of SF galaxies and AGNs. 

Key words: methods; numerical - galaxies: statistics - cosmology: theory - cosmology: 
large scale structure of universe 



INTRODUCTION 

Determining the factors that govern star formation in a 
galaxy is an important issue which is expected to shed light 
on our understanding of how galaxies are formed. It is be- 
lieved that there are several interlinked factors like turbu- 
lence, magnetic fields, cosmic rays, steller winds, supernova 
explosions, etc. functioning in a galaxy's interstellar medium 
(IS M) which together regulate star formation in the galaxy 
(see iLarsonl |2003| . I McKee fc Ostrikej |2007| for recent re- 
views). There is now mounting evidence that in addition 
to these factors operating inside a galaxy, its star forma- 
tion is also influenced by external factor s. It is quite clear 
that interactions betwee n galaxies (e.g. iBvrd fc ValtonenI 
I1990I . iMoore et al.l 119991 ) and the inte ractions of a galax; 



other galaxy proper t ies like luminosity (e.g. lEinasto et al.l 
20031. lEinasto et all l2005l. iPark et"aLl |2007| ). colour (e.j 



Hogg et al. 2OO0, Blanton et al. 2003) and morphology (e.g. 
Drossier 1980, Goto et al.|i2003i ) also exhibit environmental 
dependence. 



The analysis of filamentary patterns in the 

galaxy distribution has a long history dating back 

to a few pa pers in the lat e-seve nties and mid- 
Joeveer et al.l lll978h. 



ej.ghties 



IZel'dovich . Einasto fc ShandarinI 

Sh andarin fc ZeldovichI 1 19831 ) and lEinasto et al.l 



Einasto et al.l 



Filaments are the most striking visib l e patt e rns in the galaxy 



with the external ambient medium 



e^ 



Gunn fc Gott 



IZabludoff et~aLl 1 19961 . iPorter fc Ravchaudhurvl 120071 ) 



197: 



trigger star formation. There is also observational evidence 
that star formation is suppressed in regions where the 
galaxy density is h igh (JLewis et al.ll2002l : iGomez et al.ll2003l : 
iBalogh et al.ll2004) . This is possibly related to the fact that 



is 



distri b ution 

19961. ISha ndarin fc Ycss 



Geller fc Huchral 19891. Shectman et al.l 
Bharadwai et al. 2000 



1998, 



Miiller et a l. 
20011. 



2000, Basilakos, Plionis, fc Rowan-Ro binsop 
iDoroshkevich et all I206T 



Pimbblet. Drinkwater fc Hawkrigdl2004l ). The filamentarity 
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is found to be statistically s ignificant up to length-scales 
80fe~ ^ Mpc and not bey ond (Bharad wai. Bhavsar fc ShethI 
I2OO4I: iPandev fc Bharadwa i 2005). Our earlier work 
i Pandev fc Bharadwail I2OO6I . hereafter Paper A) shows the 
degree of filamentarity to depend on physical properties 
like the luminosity, colour and morphology of the galaxies. 
In the present work we investigate the relation between 
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log{[Nll]/iS5a6/H;j) 

Figure 1. BPT diagnostic diagram showing the ratios 
[OIII]/Hi3 and [NII]/Ha for all galaxies with Ixjai^ > 
2. The theoretical curve s eparating SF and AGNs (upper- 
dashed, lKewle^_et_alJj200l|) and the empirical curve (lower- 
dotted, iKauffmann et al.ll2003l) are shown. 



the filamentarity observed in the galaxy distribution and 
ongoing star formation activity in the galaxies. 

A brief outline of our paper follows. Section 2 de- 
scribes the data and method of analysis, our results and 
conclusions are presented in Section 3. We have used a 
ACDM cosmological model with ^mo = 0.3, JIao ~ 0.7 and 
h = 1 throughout. 



DATA AND METHOD OF ANALYSIS 



The Sloan Digital Sky Survey (SDSS. lYork et al.ll2000l ) is 
a five-passband (u,g,r,i,z) imaging and spectroscopic survey 
of the Northern Galactic hemisphere to a limiting Petrosian 
r band magnitude r < 17.77. Our analysis is limited to 
seven thin strips on the sky dr awn from the SDSS DR5 
l|Adelman-McCarthv et al.l l200p). each spanning 90° in A 
a nd 2° in T? . Here A and rj are survey co-ordinates defined 
in lStoughton et al.l l|2002l '). These strips are identical in sky- 
coverage as the ones used in Paper A and are shown in Fig- 
ure 1 of that paper. Only galaxies with extinction corrected 
Petrosian r band magnitude in the range 14.5 < rrir < 17.77 
were used. Volume limited samples are constructed in the 
same way as discussed in Paper A. The samples cover r- 
band absolute magnitude range —21 < Mr < —20 and red- 
shift range 0.043657 < z < 0.114635 which correspond to 
130/i~^Mpc to 335/i~^Mpc comoving in the radial direc- 
tion. Finally we have 17225 galaxies distributed in seven 
strips. 

Seven emission lines {Ha{6565A°), H^(4863A°), 
O/(6302A°), OJ/(3727yl°), OIII{5008A°), iVJJ(6585A°), 
SII{6718A°)) are required to classify a galaxy as either star 
forming (SF) or AGN. Only the galaxies having all these 
seven emission lines with I\/aj^ > 2 for all the lines are con- 
sidered for classification as SF or AGN, galaxies which fail to 
meet this criteria are referred to as Other galaxies. Here I\ is 
the emission line fiux and aj^ is it's uncertainty. We have fur- 
ther classified the Other galaxies as ei ther red {u — r> 2.22 ) 
or blue using the criteria proposed bv lStrateva et al.l l|200ll ). 
This color selection criteria ensures that the red galaxies 



are 'ellipticals' with 90% completeness. Most of the galaxies 
(~ 80%) classified as Other are found to be red galaxies. 

Both the SF galaxies and AGNs show strong Ha emis- 
sion. T o differentiate betwe en them we use the BPT di- 
agram IjBaldwin et al.l Il98ll ) where the logarithm of the 
ratio's of [0777]/H^ and [Nil] /Ha are plotted (Figure 
H. We used toge t her th e demarcation curve provided by 



Kauffmann eT^ (|2003h . log(^) > 0.61/(log(^) 



0.05) -|- 1.3 and the theoretical separation curve log( 



OIII\ 



0.61/( log(-|g-) - 0.47) + 1.19 provided by iKewlev et af] 
l|200ll ). Galaxies which lie below both these two curves are 
classified as SF whereas those which lie above both the 
curves are classified as AGNs. Galaxies living in the interme- 
diate region of these two curves are discarded from further 
analysis. 

Table I. gives a detailed break-up of the composition of 
the seven slices. When comparing the filamentarity of two 
different classes of galaxies it is necessary that the galaxy 
number density be the same for both the classes (Paper A). 
We have ensured this by culling the class of galaxies which 
have a larger number density. 

The Ha line was used to determine the star formation 
rate (SFR) of the S F galaxies in units of {MQyr~^) using 
the relation given bv lHopkins et al.l l|2003l ) 



SFRHaiMeyr-') = AirOfSH. 



10 



'0.4(rp 



1.27 X 10^4 



^'''""■' , Sho/Sh^ 2.114. 



2.86 



where Di is the luminosity distance, Shu and Sh/} are the 
steller absorption corrected Ha and H/s fluxes respectively 
and Vpstro and rfiter are r-band Petrosian and flber mag- 
nitudes respectively. The reader is referred to iKennicuttI 
( 1998 ) for various other definitions of the SFR. The SF 
galaxies are further classified as either high SFR or low SFR 
galaxies using the criteria SFR > 2.7 MQyr~^ which has 
been chosen so that the number density of high and low SFR 
galaxies are nearly the same. 

All the strips that we have analyzed are nearly two di- 
mensional. The strips were all collapsed along the thickness 
(the smallest dimension) to produce 2D galaxy distrib utions. 
We use the 2D "Shapefinder" statistic (IBharadw ai et al.l 
120001 ) to quantify the average filamentarity of the pat- 
terns in the resulting galaxy distribution. A detailed 
discussion is presented in Paper A , and we present 
only the s alient features h ere. The read er is referred to 
ISahni. Sat hvaprakash.fc Sha ndarinI ([1993) for a discussion 
of Shapefinders in three dimensions. 

The galaxy distribution is represented as a set of Is 
on a 2-D rectangular grid of spacing 1 /i~^Mpc x 1 h~^Mpc, 
empty cells are assigned a value 0. We identify connected 
cells with a value 1 as clusters using the 'Friends-of-Friend' 
(FOF) algorithm. The filamentarity of each cluster is quan- 
tified using the Shapefinder J^ defined as 



J^ = 



{P^ - 16S) 
{P - 4/)2 



(1) 



where P and S are respectively the perimeter and the area 
of the cluster, and I is the grid spacing. The Shapefinder T 
has values and 1 for a square and filament respectively, 
and it assumes intermediate values as a square is deformed 
to a filament. We use the average filamentarity 
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Table 1. The definition and composition of the seven strips analyzed here. 



Strip Lambda Range 



Eta Range Total Number of galaxies(All) Star forming AGN Other Other Blue 



Strip 1 


-50 < A < 40 


9 < »? < 11 


Strip 2 


-50 < A < 40 


11 < r; < 13 


Strip 3 


-60 < A < 30 


13 < r; < 15 


Strip 4 


-60 < A < 30 


15 < r; < 17 


Strip 5 


-50 < A < 40 


21.5 <») < 23.5 


Strip 6 


-50 < A < 40 


24 < »; < 26 


Strip 7 


-50 < A < 40 


26 < ») < 28 



2922 
2470 
2296 
2300 

2454 
2541 
2242 



774 


189 


1959 


333 


655 


156 


1659 


327 


649 


153 


1494 


315 


644 


133 


1523 


266 


702 


152 


1600 


283 


671 


181 


1689 


325 


617 


173 


1452 


284 



F2 






(2) 



to asses the overall filamentarity of the clusters in the galaxy 
distribution. 

The distribution of Is corresponding to the galaxies is 
sparse. Only '^ 1% of the cells contain galaxies and there 
are very few filled cells which are interconnected. As a con- 
sequence FOF fails to identify the large coherent struc- 
tures which corresponds to filaments in the galaxy distri- 
bution . We overcome this by successively coarse-graining 
the galaxy distribution. In each iteration of coarse-graining 
all the empty cells adjacent to a filled cell are assigned a 
value 1. This causes clusters to grow, first because of the 
growth of individual filled cells, and then by the merger of 
adjacent clusters as they overlap. Coherent structures ex- 
tending across progressively larger length-scales are identi- 
fied in consecutive iterations of coarse-graining. Finally a 
transition from many individual structures to an intercon- 
nected network is found to o ccur at a filling factor 0.5 — 0.6 
l|Pandev fc Bharadwaill2006l ). 

So as not to restrict our analysis to an arbitrarily chosen 
level of coarse-graining, we study the average filamentarity 
after each iteration of coarse-graining. The filling factor FF 
quantifies the fraction of cells that are filled and its value in- 
creases from ~ 0.01 and approaches 1 as the coarse-graining 
proceeds. We study the average filamentarity F2 as a func- 
tion of the filling factor FF (Figure [2]) as a quantitative 
measure of the filamentarity at different levels of coarse- 
graining. The values of FF corresponding to a particular 
level of coarse-graining shows a slight variation from strip 
to strip. In order to combine and compare the results from 
different strips, for each strip we have interpolated F2 to 7 
values of FF at an uniform spacing of 0.1 over the inter- 
val 0.05 to 0.65. Coarse-graining beyond FF ~ 0.65 washes 
away the filaments and hence we do not include this range 
for our analysis. 

Our method of analysis is similar to the Friend-of- 
Friends method with varying linking length and also the 
density field (DF) method with a fixed smoo t hing l ength 
and varying threshold density IjEinasto et al.l [2007g). All 
these methods identify density enhancements in the cos- 
mic web. With increasing lengthscale we progress from clus- 
ters (~ 1 - lOMpc) to superclusters (~ 10 - 100 Mpc) 
and finally an infinite interconnected network, the cos- 
mic web. Often individual one dimensional structures are 
identified as filaments using a variety of criteria, for ex- 
ampl e a chain of galaxies connecting two adjacent clus- 
ters IjPimbblet. Drinkwater fc Hawkrigd |2004| : IStoica et al.l 




0.6 0.2 

Filling Factor (FF) 

Figure 2. Average Filamentarity as a function of Filling Factor 
comparing the different classes of galaxies indicated in the panels. 



|2007|). Instead of focusing on the properties of individual 
structures, we quantify the overall filamentarity of the en- 
tire galaxy distribution. This is done over a range of filling 
factors. Figure 9 of Paper A shows that the average length 
of the structures identified by our method is ~ 10 Mpc at 
FF < 0.1 corresponding to galaxy clusters, < 100 Mpc for 
FF < 0.4 corresponding to superclusters and the average 
length is > 100 Mpc for 0.5 < FF < 0.6 where we have the 
percolation transition. We note that th e percolation thresh- 
old i s FF = 0.12 in the DF method (JEinasto et all |2005| . 
|2006| . l2007a ) indicating that the structures identified by our 
method are inherently thicker enclosing larger empty regions 
relative to the DF method. Despite this, it is justified to refer 
to these structures as filaments because the average filamen- 
tarity F2 of these structures is quite high (Figure [SJ for a 
large range of the filling factor FF. 



RESULTS AND CONCLUSIONS 

We first compare the average filamentarity of the star form- 
ing galaxies with that of the Other galaxies. For both classes 
of galaxies we use the results from the seven different strips 
to compute the mean and variance of the average filamen- 
tarity F2 at uniformly chosen values of the filling factor FF. 
The results are shown in the top left panel of Figure [2] We 
find that the SF galaxies have a higher F2 compared to the 
Other galaxies at all values of FF except at the smallest 
values FF < 0.1 where the Other galaxies have a larger F2. 
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At each value of FF we use the Student's t-test to de- 
termine whether the differences in F2 between the SF and 
Other galaxies is statistically significant or not. For each 
class of galaxies we have seven strips which we use to cal- 
culate the mean and variance of F2. We use the Student's 
t-test to test the null hypothesis that the mean F2 is the 
same for the two classes of galaxies. The variance in F2 
is similar for the SF and the Other galaxies, and we use 



sd 



E.e^'-'-"^)^+E. 



A^i-^Br 



(1 



+ 



to esti- 



Na+Nb-2 ^Na ' Nb- 
mate the standard error for the difference in the means. 
Here A and B refer to the SF and Other galaxies respec- 
tively, xa and xb refer to the mean F2 for the respective 
class of galaxies and Na = Nb = 7 respectively refer to 
the number of data points (different slices) for each class 
of galaxies. We use t = ^'^"'^^ to estimate the significance 
of the differences in the means. This is expected to follow 
a Student's t-distribution with 12 degrees of freedom. We 
accept the difference in the means as being statistically sig- 
nificant if the probability of its occurring by chance is less 
than 5%, Here we find that the difference in the mean av- 
erage filamentarity between the SF and Other galaxies is 
statistically significant at all values of FF. 

In an earlier work (jPandev fc Bharadwaill2006l ) we have 
studied the colour and morphology dependence of -F2. We 
find that the blue galaxies have a larger value of F2 at 
nearly the entire range of FF except at the smallest val- 
ues {FF < 0.2) where the red galaxies have a larger F2. A 
similar behaviour was also seen when the spirals were com- 
pared with the elliptical galaxies. It was found that at large 
FF the spirals have a larger F2 as compared to ellipticals, 
and the behaviour is reversed at {FF < 0.25). Most of the 
SF galaxies are blue while the Other galaxies are a mix- 
ture of blue (spiral) and red (elliptical) galaxies. As noted 
earlier, the Other galaxies are predominantly red ellipticals 
(~ 80%). To test if the observed difference in the filamentar- 
ity between the SF and Other galaxies is a consequence of 
the fact that the SF galaxies are predominantly blue galax- 
ies whereas the Other galaxies are predominantly red, we 
have separately compared the SF galaxies with blue galax- 
ies drawn from the Other galaxies. The results are shown in 
the top right panel of Figure [2] We find that the blue galax- 
ies drawn from the Other galaxies have a higher F2 than the 
SF galaxies for the entire FF range, though the difference is 
statistically significant at only the smallest FF value. The 
error-bars for the comparison of SF and blue Other galaxies 
are relatively larger than those for the comparison of SF and 
the entire sample of Other galaxies. The large error-bars do 
not permit us to rule out the possibility that the excess fil- 
amentarity of the SF galaxies relative to the Other galaxies 
is a reflection of the differences in the filamentarity of red 
and blue galaxies. 

We next compare the filamentarity of SF galaxies with 
that of AGNs. The results are shown in the bottom right 
panel of Figure [5] We find no statistically significant differ- 
ence in the filamentarity of SF galaxies and AGNs except at 
the lowest value FF = 0.05 where the AGNs have a larger 
filamentarity compared to the SF galaxies. 

Finally we divide the SF galaxies into two classes with 
equal number of galaxies based on their star formation rates 
(SFR) and compare the filamentarity of the galaxies with 
low SFR (< 2.7M0yr"^) with those with high SFR. The 
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Figure 3. Absolute magnitude distribution of high SFR (thick 
blue) and low SFR (thin black) galaxies. 



results are shown in the bottom left panel of Figure [2] We 
find that there is a statistically significant difference over 
the range 0.25 < FF < 0.45 corresponding to superclusters. 
The low SFR galaxies show a higher filamentarity compared 
to the high SFR galaxies. The filamentarity of high and low 
SFR galaxies do not show any statistically significant differ- 
ences at low and hi gh FF. 

Earlier works (IPandev fc Bharadwai1l2005l . [iooil ') show 
that the low luminosity galaxies have a higher filamentarity 
compared to the high luminosity galaxies. (Figure 7, Pa- 
per A). This is very similar to the drop in the filamentarity 
of high SFR galaxies compared to the low SFR ones. A pos- 
sible explanation for the difference in filamentarity observed 
between high SFR and low SFR galaxies is to assume that 
the SFR is related to luminosity, with the more luminous 
galaxies having a higher SFR. To test if there actually is 
a relation between the SFR and the luminosity, we have 
separately considered the luminosity (absolute magnitude) 
distribution of the galaxies in the low SFR and high SFR 
classes (Figure |3]). We find that there is a relatively larger 
number of low luminosity galaxies in the the low SFR class 
as compared to the high SFR class. Similarly there is a larger 
number of high luminosity galaxies in the high SFR class in 
com parison with the lo w SFR class. 

JBrinchmann et al.l (J2004) have used a large sample of 
star forming galaxies drawn from the SDSS to study how 
the SFR depends on various physical parameters of the in- 
dividual galaxies. In Figure 17 of their paper they show 
evidence for a strong correlation between the stellar mass 
and the SFR. Over a large range of stellar masses 6 < 
log(M«/MQ) < 10 they find that the mean SFR increases 
with M-t. The correlation breaks down at log(i\f*/M0) > 10. 
Further, they find that the strong correlation between M* 
and the SFR is a recurring theme throughout their analysis. 

Environmental effects are another important factor 
which could be responsible for the high SFR galaxies having 
a lower filamentari ty. Our ea rlier works us ing N-body sim- 
ulatio ns (Bharadw aj fc Pande v 2004; Pan dev fc Bharadwail 
[20071) show that the filamentarity depends on the bias. The 
filamentarity falls if the galaxies have a high bias ie the 
galaxies preferentially inhabit high density regions. It is thus 
possible to explain our findings if we assume that the high 
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SFR galaxies preferentially inhabit denser regions compared 
to the low SFR ones. 

In summary we have identified two possible explana- 
tions for the fact that the high SFR galax;ies have a less 
filamentary distribution than the low SFR ones: (1.) A cor- 
relation between the SFR and individual galaxy proper- 
ties like luminosity (2.) A relation between the SFR and 
environmental effects like the density. Possibly both of 
these are interconnected and are simultaneously at play. 
It is interesting to discuss these two possibilities in the 
light of other related findings. It is well known that the 
more lumin ous galaxies preferential l y inhabit denser re- 
gions (e.g. lEinasto fc Einastc Il987l . lEinasto et all |2003| . 
lEinasto et al.l 120061 . IPark et al.l 120071 .) and exhibit stronger 
clustering stren g th compared to t he fainter ones (e .g. 

TT. •^ i . ' '-ir^oo' 't~v _ ■ _ . i . 1 ' ' -1 



Hamihon' '1988", 'Davis et al." '1988", 'White, TuUy & Davifl 
1988 , Park, Vogelcy, GcUcr fc Huchra 1994, Loveday et aL 
19951 . iGuzzo et al] 1 19971 . iNorberg et all 1200 ll . IZehavi et all 
20051 ). This e ffect is quant i fied t hrough a luminosity 



bias relation llBenoist et al 



leg t nrougn a mrn i nosity 
19961: iNorberg et all I2001I : 



iTegmark et al.l |2004|: iPandev fc Bharadwail |2007|). Earher 



studies of the env i ronment dependenc e ( Lewis et al.ll2002l : 
IComez et al.ll2003l : iBalogh et al.l 12004 ) all show a suppres- 
sion of star formation activity in high density regions. While 
these works all find a decrease in the fraction of SF galax- 
ies within a distance of 2 — 3 virial radius from the cen- 
ter of clusters, th e effect of de n sity o n the SFR distribu- 
tion is not clear. iGomez et aD (J2003l ) find that the SFR 
distribution is strongly shifted to lower values in high den- 
sity environments, the effect being mos t pronounced f or the 
strongest star forming galaxies whereas iBalogh et al.l (|2004l ) 
find t hat the SFR d i stribut ion is independent of environ- 
ment. lEinasto et al.l l|2007ah find that the fraction of SF 
galaxies in rich groups/clusters in high density regions of 
rich supercluster is smaller than what is found in poor su- 
perclusters in the field. In their analysis superclusters are 
connected non-percolating systems with densities above a 
certain density threshold. lEinasto et al.l l|2007lj ) also find 
that the more luminous and richer superclusters have a 
higher degree of fil amentarity compared to the poor ones. 
lEinasto! (Il99ll ) and lEinasto et "all (|2007al ) show that both 
the local and the global environment are possibly impor- 
tant in infiuencin g morpho l ogical properties and star for- 
mation. Recently iLee fc Lil (|2008l ) have studied the corre- 
lation between the large scale environmnet of galaxies and 
their physical properties in the SDSS and the 2Mass Red- 
shift Survey. They find that the physical parameters re- 
lated to the recent star formation history are linked to the 
shear of the large scale environment of galaxies. A study 
of star formation along the Pisces-C etus Supercluster fila- 
ments (|Porter fc Ravchaudhurvll2007l ) finds that though the 
SFR and the fraction of SF galaxies declines steadily to- 
wards the cores of clusters, there is an increased SF activity 
in a narrow distance range around 1.5 to 2 times the virial 
radius of the cluster involved. A study of the environment 
and c lustering propertie s of star-burst galaxies in the 2dF- 
GRS (lOwers et al.l 120071 ) finds that a significant fraction of 
star-burst galaxies show morphological evidence for ongoing 
or recent tidal interaction or merger. 

There appears to be no consensus on how star forma- 
tion depends on the environment and some of the findings 
appears to be at odds with our results, we note that the as- 



sociated length-scales are quite different. While the earlier 
studies probed groups and clusters of galaxies ie. length- 
scales less than a few Mpc, the structures identified by our 
method have lengths of the order of 100 Mpc/h and thick- 
ness of the order 5 Mpc/h (Figure 9, Paper A) for the range 
of FF where the difference in filamentarity is statistically 
significant. Our interpretation is consistent with the recent 
analysis bv lShiova et al.l (|2007l ) who has shown that the high 
SF galaxies at z « 0.24 are more strongly clustered than the 
low SF galaxies. 
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